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Abstract 

The Nordstrom- Vlasov system describes the dynamics of a self-gravit- 
ating ensemble of coUisionless particles in the framework of the Nordstrom 
scalar theory of gravitation. We prove existence and uniqueness of clas- 
sical solutions of the Cauchy problem in three dimensions and establish 
a condition which guarantees that the solution is global in time. More- 
over, we show that if one changes the sign of the source term in the field 
equation, which heuristically corresponds to the case of a repulsive grav- 
itational force, then solutions blow up in finite time for a large class of 
initial data. Finally, we prove global existence of classical solutions for 
the one dimensional version of the system with the correct sign in the field 
equation. 



1 Introduction 

In astrophysics, systems such as galaxies or globular clusters are often modeled 
as a large ensemble of particles (stars) which interact only by the gravitational 
field which they create collectively. In such systems collisions among the par- 
ticles are sufficiently rare to be neglected. Let f = f{t,x,p)>0 denote the den- 
sity of the particles in phase-space, where tsM denotes time, xSR'^ position, 
and p G M'^ momentum. This density function satisfies the Vlasov equation — a 
continuity equation on phase space — coupled to the field equation(s) for the 
gravitational field. In the non-relativistic case this results in the Vlasov-Poisson 
system, in the relativistic setting one obtains the Einstein- Vlasov system. While 
for the former system the initial value problem is by now well understood — cf. 
[3 El ^1 EI — this is not so for the much more difficult Einstein- Vlasov system, 
the difficulty arising from the different nature of the field equations. We refer to 
^ for a recent review on the latter system, cf. also |11[I12) . In the present paper 
we investigate a different relativistic model which is obtained by coupling the 
Vlasov equation to the Nordstrom scalar gravitation theory In this theory. 
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the gravitational effects are mediated by a scalar field cj), and the system reads 

,40 f, dp 



d 



'A, 



5tf + p-V:,f- 



(at0+p-v.</))p+(i+p2)-i/2v,</) 



•Vpf = 0. 



(1.1) 



(1.2) 



Here 



denotes the relativistic velocity of a particle with momentum = and 
units are chosen such that the mass of each particle, the gravitational constant 
and the speed of light are all equal to unity. A solution (f , (f) of this system is in- 
terpreted as follows: The space-time is a Lorentzian manifold with a conformally 
flat metric which, in the coordinates (t,x), takes the form 

<7p, = e2*diag(-l, 1,1,1). 

The particle distribution defined on the mass shell in this metric is given by 



fph(i,a;,p) = f(i,2:,e'^p). 



(1.3) 



More details on the derivation of this model are given in '2'. It should be 
emphasized that the system does not constitute a physically correct model. 
Our interest in it is motivated by the fact that while being much simpler than 
the "correct" relativistic model, the Einstein- Vlasov system, it still captures 
some of the essential features which distinguish the latter from non-relativistic 
models. In particular, the Nordstrom- Vlasov model does allow for propagation 
of gravitational waves, even in the spherically symmetric case. The hope is that 
the analysis of this model will lead to a better mathematical understanding 
of a whole class of nonlinear partial differential equations and eventually to a 
better understanding of the Einstein- Vlasov system. However, scalar fields do 
play a major role in modern theories of (classical and quantum) gravity, cf. 
|3] and references therein, as well as in numerical investigations of relativistic 
gravitational waves, ^1] . 

It turns out to be convenient to rewrite the system in terms of the new 
unknowns (/,<?ii), where / is given by 



/(t,x,p) = e4*(*'-)f(t,x,p). 
The system then takes the form 



Kt^x)^ / fit,x,p) 



dp 



Sf- [(50)p+(l+/)-i/2v,J •Vp/ = 4/50, 



(1.4) 

(1.5) 
(1.6) 

(1.7) 
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where 



s = dt+p-s/: 



is the free-transport operator. The function fj, is the trace of the energy mo- 
mentum tensor. In the following we refer to the system (|1.5|) . 1)1. 6|l . (|1.7|l as the 
Nordstrom- Vlasov system. We supply it with the initial conditions 

f{0,x,p) = r{x,p), (/.(O,x) = 0i,"(x), 9t(/)(O,x) = 0i"(x), x,peR\ (1.8) 

A system of equations which belongs to the same class of kinetic models and 
is used in plasma physics is the relativistic Vlasov-Maxwell system. For this 
system Glassey and Strauss proved that classical solutions to the initial value 
problem with regular data exist and no shocks or singularities form as long as 
the maximal momentum of the particles is controlled, cf. 4 . One purpose of 
the present paper is to establish an analogous condition for global existence of 
classical solutions for the Nordstrom- Vlasov system. In order to state the result 
in a precise way, let us consider the quantity 

r{t) = sup{\p\:0<s<t, (a;,p)esupp/(s)}, 

the maximum momentum of the particles up to the time t; f{s) will always be 
compactly supported. In addition, we need to monitor the quantity 

Q(<) sup{|0(s,2;)| : < s < i, G supp/(s)}, 

and we define A{t) — 1^(1) + Q{t); for t = we define these quantities by what is 
induced by the initial data. 

Theorem 1 Initial data f'^eC^iR^), 0j,"eCf(M3), (jji^ (z (jgs-) f'"'>0 
launch a unique classical solution (/,<?!)) G C^([0,Tniax[xR®) x C^([0,Tmax[xM^) 
of the Nordstrdm-Vlasov system i) j.5|) . i) j.6|) . ji. 7| ) in a maximal interval of 
time, satisfying the initial conditions fl.8\) . // A(Tmax) < oo then T'max = oo, 
i.e., the solution (/,(/)) is global in time. 

We restrict ourselves to study the Cauchy problem forward in time; the system 
is time-reversible. The notation on the function spaces above is standard: The 
index c refers to compactly supported functions, and the index b means that all 
derivatives up to the indicated order are bounded. This boundedness condition 
which only aims to make our estimates more transparent can be removed using 
the finite propagation speed which in turn follows from local conservation of 
energy, see below. We observe that unlike the case of Q where global existence 
follows from a bound on V{t) alone we require a bound on in the support 
of / in addition. This will imply a bound on /; as opposed to the Vlasov- 
Maxwell system this bound on / is a priori not obvious since / is not constant 
along characteristics of the Vlasov equation. Note however that by (|1.3|l the 
"physical" momentum of the particles is obtained by rescaling p by the factor 
e'^ so that the meaning of Theorem ^ remains the same as with the analogous 
result for the Vlasov-Maxwell system: The occurrence of singularities in the 
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Nordstrom- Vlasov system — if any — can only be due to particles which pick up 
arbitrarily large momentum as the system evolves. 

The proof of Theorem ^ is presented in Sections |3HS1 ^ key ingredient being 
certain representation formulas for the first and second order derivatives of </> 
which we give in Section 0] These formulas are similar to the ones introduced 
by Glassey and Strauss for the relativistic Vlasov-Maxwell system, and we will 
be able to reduce many details to the case of 0. 

An interesting mathematical feature of the Nordstrom- Vlasov system is that 
if we change the sign in the right hand side of (|1.5|l then solutions blow up in 
finite time for a large class of initial data: 

Theorem 2 Consider the system obtained by coupling i) j.6|) . j j. 7| ) to the equa- 
tion 

instead of lil.5\) . Then the assertion of Theorem^ remains valid. If the initial 
data are homogeneous (i.e. independent ofx) in a ball Bfj(0) — {x (EM.^ : \x\ <R} 
for some R>0 and if 



i?V(0,0)>2, 0i"(O)>V2MO;O), 



This result will be shown in Section |3 It is interesting because the blow-up 
results from the nonlinearity in the original field equation, which has only been 
hidden by rewriting the system in the form (|1.5|l . 11.6|l . (|1.7(l . and which is 
a feature that distinguishes this system from the relativistic Vlasov-Maxwell 
system where the field equations by themselves are linear. 

Solutions of the Nordstrom- Vlasov system preserve energy, which is positive 
definite in the situation of Theorem ^ but indefinite in the situation of The- 
orem 121 For the three dimensional system considered so far we are not going 
make use of conservation of energy. However, in the last section we will show 
that conservation of energy can be used to prove global existence of classical 
solutions to the one dimensional version of the system: 

d?cb-dlc^=-^i, (1.9) 

fi{t,x)^ f f{t,x,p)^^, (1.10) 

Sf- [(50)p+(l+p2)-i/2a^^] dpf^2fS<l>, (1.11) 

where now x,pGM.. The right hand side of the Vlasov equation modified so as 
to retain conservation of energy, cf. the next section. Again, the sign in the field 
equation turns out to be important. While the blow-up result of Theorem |21 
also holds in other dimensions we have the following global existence result for 
the "correct" sign: 

Theorems Initial data f"'eC^{R^), e (R) , <j)^'eCl{R) with f">0 
launch a unique global classical solution (/,(/>) G C"'"([0,oo[xM^) x C^([0,oo[xM) 
of the system ^TTBil, 
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2 General properties of the Nordstrom-Vlasov 
system in spatial dimensions 

In this section we collect some general properties of the system under considera- 
tion. Some of these will be important in the sequel while others are included for 
reasons of general interest. Since in Sectional we consider the one dimensional 
version of the Nordstrom-Vlasov system, we discuss here the generalization to 
the case of N spatial dimensions, N>1. This generalization is defined as the 
system obtained by adapting the derivation of |2] to the case in which the space- 
time is a 1 + N dimensional Lorentzian manifold. This assumption affects the 
argument of 2 in a non-trivial way only in two points. Firstly, the modulus 
of the determinant of the physical metric, denoted by g in 0, changes into 
g = exp[2(l4-iV)(/))]. Secondly, the rescaling law which relates the stress-energy 
tensor in the physical and unphysical frames becomes, using the notation of 
12], Tr = e-Kp[{N + 3)(l)]Tf'. The latter identity follows by the definition of the 
stress-energy tensor, which is independent of the frame and of the dimension, 
as the variation with respect to the metric of the matter action. With these 
modifications one can show that the generalization of (|1.5|) , (|l.t)|l , H1.7|l in 1-1- A*" 
dimensions is obtained by taking i e M, xE R^, p G and replacing the Vlasov 
equation p. 7(1 by 



5/- 



(5</))p+(l+p2)-i/2v^^ .Vp/=(l + 7V)/50. (2.1) 



The resulting system reduces to (|1.5|) . (|1.6|l . (|1.7|l in the case iV = 3 and to 1(1. 9() . 
l(TTn|l . ((rTTl) for iV=l. 

An important tool in the sequel are the characteristics of the Vlasov equa- 
tion. Let {X,P){s) = {X,P){s,t,x,p) denote the solution of 



(2.2) 



dx p 
ds ' ^l+p2' 

— = -{S(j)){s,x,p)p , 2.3) 

ds VI +P 



with initial condition 



{X,P)it,t,x,p)^{x,p). (2.4) 
Integrating ((2.1(1 along characteristics we obtain 

f{t,x,p)^r{X{0),P{0)) + {l + N) f {fS<t>){T,X{T),P{T))dT. (2.5) 

Moreover ((2.11) can be rewritten as 

5(e-(i+^)^/)_ [(50)p+(l+p2)-i/2v,0] .V,(e-(i+^)^/) = 0, (2.6) 

showing that the function e"'^"*"^'"^/ is constant along solutions of 1(2.2(1 . ((2.3(1 
(which of course are the characteristics of (12.6(1 as well) , and hence the function 
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/ can be also represented as follows: 



f{t,x,p) = r{X{0),P{0))exp[{l + N)cf>{t,x)~{l + NW'\X{0))] . (2.7) 

By (|2.5|) . the function C(s) — f{s,X{s),P{s)) solves the differential equation = 
{1 + N)r]{s)(, with ri{s) — (j){s,X{s)). Integration of this equation yields 1)2. 7|) . 
The latter equation shows that the compact support property of f{t) which is 
assumed initially propagates. 

It should be noted that the characteristic flow is not measure preserving on 
M.'^^ . Indeed, the fact that the (a;.p)-divergence of the right hand side of H2.2|) . 
H2.3)l equals —NS(f> implies that 

det^^^{0,t,x,p)^exp[N<l>{t,x)-Nc^'^X{0,t,x,p))], 

which together with H2.7|l implies conservation of rest-mass: 

/(t,a;,p)e"'^(*'"=) dxdp = const. 



//■ 



The Nordstrom- Vlasov system also satisfies an energy conservation law. In the 
formulation of the system used here, the energy density is defined by 

eit,x)= J ^/l+^f{t,x,p)dp+^{dt^f + ^{Vx4>f, 

and the momentum density by 

)p{t,x)= j pf{t,x,p)dp-dt4>Vx(l)- 
It is easy to check that local energy conservation holds, 

which upon integration leads to conservation of the total energy: 



^Jl+p^f{t,x,p)dpdx+-j [(dt(t>f + {Vxcj)f]{t,x)dx = const. 

We note that the definition of energy and rest mass is independent of the dimen- 
sion. Moreover, the total energy is positive definite. However, the sign in front 
of the field energy is reversed for the "repulsive" case considered in Theorem El 
and so the total energy is no longer positive definite in that case. 

The contents of this section will be applied for = 1 in Section [7| and for 
iV = 3 in the rest of the paper without further comment. 
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3 Preliminary estimates 

We fix some notation. The usual L°°-norm is denoted by || • ||oo, || ■ lU.cxi is the 
L°°-norm of the argument and its derivatives up to order k. Moreover, 



\\Df{t)\\^ ^ sup{\d,J{t,x,p)\,\dpJ{t,x,p)\:{x,p)eM.', i = l,2,3}, 
WD^moo = sup{|9t(/.(t,x)|, \d,Mt,x)\-x€R^ 1 = 1,2,3}, 
\\D^<l>{t)\\oo = sup{\d^<j,{t,x)\, \dtd^^(l){t,x)\, \d^^d^^(l){t,x)\:xeR\ z,j = l,2,3}. 



We shall use similar definitions for the one dimensional case in Section |7| Nu- 
merical constants are denoted by c and may change from line to line. Constants 
depending on the initial data are denoted by C,Cq,Ci,..., dependence on the 
length of some time interval is denoted by Ct- Since has compact support 



without loss of generality we may assume that Pq>1 and hence 'P{t)>l for all 
t, otherwise we add 1 to both quantities. 

We prove a few preliminary estimates which will be used in the proof of 
Theorem n By solution we always mean a solution of the Nordstrom- Vlasov 
system (|1.5f) . (|1.6|l . (|1.7|l with the regularity and initial data as specified in 
Theorem ^ 

Lemma 1 For any solution defined on some interval [0,r[, 



Po = sup{|p| ■.ix,p)e supp/'"} < oo; 




// Q{T) <oo then 

sup{||/(<)||oo:0<i<T}<oo. 
Proof: The estimates are immediate by (|2.5|l and (|2.7|l . 



□ 



Lemma 2 For any solution defined on some interval [0,r[, 




where Co^c{U\^%^oo + Ut\\oo). //A(r)<oo then 



sup{|j0(O||oo:O<i<T}<oo, 



i.e., (j) is bounded everywhere and not just on the x-support of f . 



Proof: By Kirchhoff's formula, 




(3.1) 
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The first term is the solution of the homogeneous wave equation with data 
(0q",0'") and is estimated by 

||<^o(i)||oo<c(l+t)(||C||l,oo + ||<^rl|oo)-Co(l+i). 

The second term is estimated by 



'0 J \x~y\—t- 
nt 



/ fi'r,y,p) 

-tJ\p\<V{t) 



dp dr 



't-T 
ft 



<c it-r)\\fiT)\\^rirrdr<cil + t) ||/(r)||oo7'(r)"dr. 







□ 



Together with Lemma Q] this completes the proof. 
Lemma 3 The function V satisfies the estimate 

Vit)<Po + 2f \\D^{T)\\^{l+V{T))dT. 

Jo 

Proof: This follows by integrating (|2.3|) and observing that the right hand side 

-{dt(j>iT,x) +p- V,0(t,x))p- (1 +p2)-i/2v^0(^^2:) 

of that equation is estimated by 2||D(/)(r)||oo(l + b|)- ^ 
To conclude this section we estimate the derivatives of the characteristics 
with respect to their initial data in terms of the derivatives of (/). 

Lemma 4 Assume that V {T) < oo . Then for all te[0,T[ and x,peM.^ , 



|V(,,p)(X,P)(0,i,a:,p)|<CTexp 



C7 



(l + |lD0(r)|U + |li?20(^)|U)dT 



Proof: Differentiating H2.2(l and (|2.3|l with respect to Xi and integrating we 
obtain 



dpk 
t dpj 



(P(T))9,,Pfe(r)dT, 



{Sd,,(b)pj + il+p^)-'/^d.,,d,^(b\ {T,X{T),P{T))d^Mr)dT 
Pkdx,4'' 



{S(j))5jk + j—pidxi(t)-\ , 2 

Opj 1+J3^ 



{T,X{T),P{T))dx,Pu{T)dT, 



summation over repeated indexes being understood. Now 

dpk _ {^+p'^)5jk~'PjPk 



dp. 



(l+p2)3/2 
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By assumption, P{t) < Ct, and so we get, for s<t, 
|V,X(s)| <c+ / |V,P(r)|dr, 

J s 

|V,P(s)| <Ct f (||i?0(r)|loo + P'0(r)|U)(|V,X(T)| + |V,P(T)|)dT. 



If we add these estimates and use Gronwall's inequality we obtain the assertion 
for the a:-derivatives, and the argument for the p-derivatives is analogous. □ 



4 Representation and estimates for Dcj) and D'^cj) 

Let us first recall the structure of the integral representation for the electromag- 
netic field {E,B) in the case of the relativistic Vlasov-Maxwell system, cf. 01 
Thm. 3]: 

Eit,x)^EDit,x)- [ j a''{oj,p)f{t-\x-yly,p)dp- '^^ 



-y\<t 



x — y\ <.t J 



\x-y\ 



2 



~a^{uj,p){Sf){t-\x-yly,p)dp-^ (4.1) 

F y\ 



B{t,x) = BDit,x) - [ [ a^{Lu,p)f{t 

J\x^y\<tJ 



\x-y\,y,p)dp- 



\x-y\ 



2 



-y\<t 



h^{u;,p){Sm- \x-yly,p)dp-^ (4.2) 

\x-y\ 



Here uj = {y~x)/\x — y\, Ed^Bjj are integrals of the initial data and the kernels 
are given by 



(l+p2)(l+a;.p)2' ^ '^^ l + uj-p' 
(w,p) = 7-— ^r-— a (uj,p) = - 



(l+p2)(l+^.^)2' V '^^ l+LU-p 

In the next two propositions we show that the derivatives of (/) satisfy very 
similar integral representations. 

Proposition 1 The time derivative of (p satisfies the representation formula 
dt^{t,x) = dtMt,x)-t-' f f J^ dSy 



x-y\=t 



\x-y\<t 



a"^* i^^:P)fit-\x-y\,y,p)dp 

\x-yr 
dy 



b"*" iuj,p)fiS(f)){t- \x-y\,y,p)dp 



\x-y\<fJ \x~y\ 

dy 



\^-y\<t 



c"^* {uj,p)-Va;(f)fit-\x-y\,y,p)dp 



\x-y\ ' 
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where the kernels are given by 



v/l+p2(l+cj-p)2 



(l+p2)3/2(l+^.^)2- 

They satisfy the estimates 

|a'^'(u;,p)|<c(l+p2), |6^'(c^,p)|<cVT+p2, |c'^'(u;,p)| <c. 

Proof: The claim is proved by a calculation similar to the one in |4| Thm. 3]. 
We give only some of the details. By (|3.1|l . 

dt4> = Mo-t-' I I f'"\y,p)^$=dSy 



x-y\=t 



dtf{t- \x-y\,y,p) 



dp dy 



l\x-v\<tJ ^/l+p^\x-y\ 

In the last integral we use the identity 

{dtg){t-\x-y\,y,p) = {l+uj-py^ {Sg){t~\x-y\,y,p) 

-p-Vy[git~\x-yly,p)]\, (4.3) 

which holds for every function g — g{t,x,p), and integrate by parts in y to 
obtain 



dt(j) = dt<i)o{t,x)-t / / -— ^ 



-.dSy 



\x-y\<t 



x~y\<t 



\x^v\=tJ y/l+p^ 

a''"{^,p)fit~\x-y\,y,p)dp 

\x-yr 

d^^icj, p){Sf){t-\x-y\,y,p) dpj^^ , (4.4) 



where a'^' = {1+p'^)^^ {l+uj -p}^^ . Hence we see that the function has 
been written in the form ()4.1|l with replaced by a"^' and by a"^'. To 
complete the proof we observe that by the Vlasov equation Ijl.TH . Sf — F{t,x,p) ■ 
Vpf + 4:fS(j) where F{t,x,p) = {S(j))p+{l+p'^y^/^\7x(j> is the force term in the 
Vlasov equation. Substituting into the last integral of H4.4|l and integrating by 
parts in p the representation formula for dtfj) is obtained after a straightforward 
calculation. For the bounds on the kernels we use the estimates 

^ + ^-P^^^^' l^+i^'<2(l + ^-F), (4.5) 
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cf. Lemma 1]. Thus for example, 

\c^^{u,p)\<cil+p')-^l^{l + u-p)-^'^<c{lW)-'^^ 
estimates for the other kernels being similar. □ 



Proposition 2 The spatial derivatives of cf) satisfy the representation formula 

J\x-y\=tJ l+^-P Vl+P 



\x^y\<t 



\x-y\<t 



\x-y\<t 



af^ {uj,p)f{t~ \x-y\,y,p)dp '^^ 

dy 



bf-iLU,p){fS^)it~\x~yly,p)dp 



Ci iu;,p)-V^(j)f{t-\x-y\,y,p)dp 



\x~y\ 
dy 



\x-y\' 



where the kernels are given by 



a-^-^ Vl+P^[a^-pAa^], cf^ ^Wic"**, i = 1,2,3. 

They satisfy the estimates 



\a^^{oj,p)\<c{l+p^), \b^^{uj,p)\<c^/lTp', \cf^{u;,p)\<c. 
Proof: The only difference to the proof of Proposition^is that now the identity 
idy,g){t-\x-y\,y,p) = uj,{l + uj-py^{Sg){t-\x-y\,y,p) 

+ (^S,, - -^^^ dy, [git-\x-y\,y,p)] (4.6) 

is used instead of H4.3|l . Thus the analogue of 14.41) in this case is 

VMt,x) = (V.0)d- / / a^''iuj,p)f{t-\x-y\,y,p)dpj-^ 



\x-y\<t 



x — y\<t J 

~a'*'^{u;,p){Sf){t-\x-y\,y,p)dp 



dy 



\x-y\' 



where {Vxfj^jD is the data term and a'^'' = {l+p^) + ^uj. Since 

and are bounded by cy^l+p"^ (cf. the estimates on the kernels follow. 

□ 

We use the above representation formulas to estimate the first order deriva- 
tives of 0. 
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Lemma 5 For any solution on some time interval [0,T[ 



I! W)l|oo <(! + <) 



Cl+C / {l + \\D4>{T)\\^)\\J{T)\\ooV{TfdT 



ie[o,r[. 



«;;iere Ci=cP4(l + ||f"||oo + ||C||2,oo + ||</'r||i,oo). IfA{T)<^ then 

sup{||i:)0(i)||oo ■.0<t<T}<oo. 

Proof: By Proposition ^ dt4> is the sum of five terms wliicli we denote by 
/i , . . . , /s . The first term is estimated as 

|/l|<c(||Cll2,oo + ||0rl|l,oo)(l+i). 

For the second term we use (|4.5ll to obtain 



|/2i<c(i+t)iirii 



/ Vi+^dp<cP^\\n\^{i+ty, 

J\p\<Po 



recall that Po > 1- For we have, using the bounds on the kernels in Proposi- 
tion □ 

I/3I < c fvirf [ [ f{r,y,p)dpdSy-^ 



<c ||/(r)||ooP(r)'^dT. 
Jo 

Likewise we obtain 



|/4|<c(l+i) / \\D^it)\\oo\\fiT)\\ooViTfdT, 

Jo 

\h\ < C{l+t) f \\Dcl>{t)\\^\\f{T)\\^V(TfdT. 

Jo 

To complete the proof we add these estimates and repeat the same argument 
for Va;0. The bound on D(j) follows from Gronwall's inequality together with 
Lemma ^ □ 

If we differentiate H2.7|l and use the estimates from Lemmas 0] and we 
immediately obtain: 

Corollary 1 Assume that A{T)< 00. Then 

P/WI|oo<CTexp[cT / (l + p20(r)||^)dTl, <e[0,T[. 



Next we derive the integral representations and the necessary estimates for 
the second order derivatives of 6. 
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Proposition 3 The function dfcf) can be represented by the formula 



dUit,x) = {dU)D+j> I a'^"{LO,p)fit-\x-y\,y,p)dpj-^ 



x-y\<t 



-y\<t 



b'*'"{LO,p){Sf){t- \x-y\,y,p)dp 



dy 



\x-y\^ 
dy 



c^"{iu,p){S^f){t-\x-y\,y,p)dp- 

x-y\<tJ \^~y\ 



where {df4>)D contains only integrals of derivatives of the initial data and the 
kernels are smooth. Moreover, the kernel a"^" of the most singular integral 
{df4))sing has zero uj-average: 



a'^"{uj,p)dS^^Q, pel 



)3. 



observe the principal value of the spatial integral in (9^0)sing- Similar represen- 
tations with kernels having the same properties hold for the other second order 
derivatives of (j). In particular, in each case the kernel of the most singular 
integral has zero Lo-average. 



Sketch of the proof: The representation for dfcj) is obtained from (|4.4|) in the 
same way as Glassey and Strauss derive the representation for the gradient 
of the electromagnetic field from H4.1|I " H4.2|I cf. 4, Thms. 4,5]. Each term in 
H4.4|l is differentiated with respect to t and then H4.3|) is used with g = f and 
g = Sf. Integration by parts leads to the formula for dfcj) after a straightforward 
calculation. The proof of the averaging property for a*^" can be reduced to the 
case of T using the relation between the kernels for the Nordstrom- Vlasov 
system and those for the relativistic Vlasov-Maxwell system which we noted in 
Propositions Hand 121 The crucial point is that to get from and to a'^* or 
a'f'^ we need to multiply by terms depending only on p. For example, in the case 
of d^(f> the most singular integral is obtained by differentiating in time the third 
term in the right hand side of (|4.4|) . substituting (|4.3|) and then integrating by 
parts in y. Hence the resulting kernel is 



a'^"(cj,p) = -\x-y\^p-W 



a^*iuj,p) 



{l + uj-p)\x-y\'^ 



1,3 



dyj 



{l + uj-p)\x~y\'^ 
{l + LO ■pY\x — y\'^ 



= \/l+p'^p-a{uj,p), 
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where a{u},p) is the kernel which appears in (44) of 0. Thus the averaging 
property for a*^" follows from (43) of The analogous argument applies to 
the other second order derivatives. □ 

Lemma 6 Under the assumption A{T)<oo, a classical solution on the time 
interval [0,T[ satisfies the estimate 



\D'(t>moo<CT l + ln* sup ||D/(t)|| 

L 0<T<t 



\D^cl>{T)\\^dT , te[o,T[, 



where 



In* s : 



s for s < 1 
1 + lns for s> 1. 



Combining this with Corollary^ and using GronwalVs inequality implies 

sup{|| D^<?!)(0 II oo : < t < T} < c». 

Sketch of the proof: The proof is similar to 01 Lemma 6] and is outlined only for 
d^cj). The data term in Proposition |21 is estimated in terms of the appropriate 
norms of the initial data. The second term, i.e., the most singular integral, is 
estimated using the averaging property of its kernel as follows, cf. also (59)-(63) 
of SI: 



\{dU)su.^{t,x)\<CT 



t-d 



a*"(c.,p) 



{t~T)-^\\f(T)\\^dT 

f{t-T,X + TUJ,p)-f{t~T,X,p) 



dpdS^dr 



<CTHt/d) + CT sup \\Df{T)\\ood 

0<T<t 

which holds for any <d<t <T. Let ct = supQ<^<( ||Z?/(r) 
d = t, else choose d = inm{a^^ ,t}. Then 



If (T < 1 choose 



\{dUUng{t,x)\<CT 



1 + ln* sup ||D/(r)|| 

0<T<t 



For {df(f>)s we express Sf by t\1.7\\ and integrate by parts in p. By the assumed 
bound on V the resulting kernels are bounded on the support of /. Since 
the resulting expression contains only terms in / and first order derivatives of 
(j) which are bounded according to Lemmas ^ and [S] we obtain the bound on 
|(9|(/))5|. To estimate {df(l))s2 we expand S^f using the Vlasov equation: 



S^f = 45[/50] + 5[((50)p+(l+p2)-i/2v^^).v^j] 

= 4/^2^ + 4(50) (Sf) + (520)p. Vp/ + (1 +p2)- i/2(5v, 

+ F{t,x,p)Vp{Sf)-Fk^d,J. 

dpj 



14 



After substituting this expression into {d^<j>)s2 we remove the p-gradients by 
integration by parts. Then we use H4.6|l to express dx^f and expand Sf using 
the Vlasov equation. We end up with an expression which contains only first 
order p-derivatives of /, second order derivatives of 4> and the perfect derivative 
Vy[/(t— \x — y\,y,p)]. Integrating the corresponding terms by parts in p and y 
and using Lemmas and [S] again we get the estimate 



\{dU)s<t,x)\<CT 



1- 



, dr 



Summing up the various estimates and using the same argument for all second 
order derivatives proves the first estimate of the lemma. The second one then 
follows via Corollary n] and Gronwall's inequality. □ 



5 Proof of Theorem [T] 

Let us start by proving that the initial data H1.8|l with the regularity specified 
in the theorem launch a unique solution in some local time interval. For this 
purpose, we define an approximation sequence (/*•"', (/''■"^) in the following way: 
For n = we set f^°^t,x,p) = f''{x,p) and (j)^°'> {t,x) = (l)\^{x). If 
is already defined, we define the n"^ iterate as the solution of the system 

- (1 + p2)-l/2 V^0(n-1)] . ^^f(n) ^ 4j(n) ^^("-l) ^ 

with initial data 

/(")(0,x,p) = f"(x,p), 0(")(O,x) = C(x), 5t0(")(O,a;) = 0i"(a;). 
We also define 

{t) = sup I IpI : < s < i, {x,p) e supp/(") (s) } 

and fn,m — /'•"'' — Z'-™-' , 4>n,m — ~ 0'-™'' • Local existence of a classical solution 
will follow from 

Proposition 4 For all neN, f'^'"^ eC^{[0,oo[xR^), e C2([0,cx)[x]R3). 

There exists a time > and a continuous function z : [0,T»[^ [0,oo[ such that 

7,(n)(i) + ||/(n)(i)||^ + ||^(«)(t)||^ + ||i^0(n)(i)||^<^(i)^ te[0,T4, n G N. 

Moreover 

\\fn,mit)\\oo + \\4'n,mit)\\oo + \\D<pn^jnit)\\ao + \\Dfn^mit)\\oo + \\D'^^^^ 

as n^m~f oo, uniformly on any compact time interval [0,r] C [0,r,[. 
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Proof: We adapt the estimates proved above for solutions of the Nordstrom- 
Vlasov system to the approximation sequence [f^^\(jP^'^). First we notice that 
the functions ^t^*-"-* and Vx^*-"-* have integral representations analogous 

to the ones proved for the solution: 

</)(") (i,x) = </>o(t, a:) - / [f^-\t-\x-y\,y,p) (5.1) 

J\x-y\<tJ ^/l+p^\x-y\ 



J\x-y\=tJ (1+W-p) ^l+p2 



x—y \ <.t 
\x-y\<t 



a^^{u:,p)f^^\t~\x-yly,p)dp '^^ 



\x-y\^ 



6^*(a;,p)[(50("-i))/(")](i-|x-2/U,p)dp 
c^* {u,p) ■ V,0("-i) /(") {t-\x~yl y,p) dp 



dy 



\x-y\ 

dy 



j\x~y\<tJ \x~y\ 

and similarly for Vj;(/)*^"K Moreover, if we denote by (X„,P„) the solution of 
(Ei2l-lE3|l with (j> replaced by then we have, as in (|2.5|l and (|2.7|) . 

/(")(i,a;,p) = r(X„(0),P„(0)) + 4 f\f^^^ S^'^'^-'^]{TMT),P.^{T))dT, 

Jo 

f^^\t,x,p) = r(X„(O),P„(O))exp{40("-i)(t,a;)-40'"(X„(O))}. 
Thus the analogues of the estimates in Lemmas |21 and |5l read 

||0("Hi)l|oo < Co{l+t)+c{l+t) f\\f^''\T)\\ooV^-\TfdT, 

Jo 

P0("ni)lloo <Ci(l + i) 

+ c{l+t) Al + ||i?0("-i)(T)||o,)|l/(")(r)|U7'(")(r)5dT 
Jo 

and hold for all t>0. The analogue of the first estimate in Lemma ^ is 

||/^"'Wlloo<||/"||oo+4 /*||i?^("-i)(r)|U||/(")(r)|Udr. 
Jo 

Finally, as in Lemma |31 

r^-^^t) <Po + 2 f ||i?</)("-i)(r)|U(l+7'<"Hr))dr. 



Let us define 

Qn{t) =V^"^ (t) + II /(") (t) lU + ||0(") (i) II oo + II (i) II 
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Combining the previous estimates we obtain the incquahty 

Qn{t)<C + C f Qn-l{T)Qn{Tf dT, ^^[0,1], 

where 

C = nmx{4,2Co,2Ci}p4(i + Po + ||/in||^ + ||0j,n||2_^ + ||0in||,_^). (5.2) 
It is then straightforward to show by induction that 

Qn{t)<z{t), te[0,T4, neN, 

where z is the maximal solution of 

i^Cz\ z(0) = C, 

which exists up to time r* = C^^/6<l. Since T''-"' is bounded, the regularity 
of the approximation sequence follows by (|5.1() and the induction hypothesis. 

In any time interval [0,r]c[0,r*[ and for any iterate, the quantity A(r) 
is bounded by a constant Ct- By Corollary ^ and Lemma El applied to the 
iterates we obtain uniform bounds on the first order derivatives of /^"^ and the 
second order derivatives of 0^"^ on [0,T]. The Cauchy property for and for 
</)("^ and its first order derivatives is then easily deduced by deriving a Gronwall 
inequality for fn,m- The Cauchy property for the next level of derivatives is 
less simple, and we refer to the analogous arguments for the relativistic Vlasov- 
Maxwell system, cf. gl pp. 82-85]. □ 

To prove uniqueness of local solutions let {fi,4>i) and (/2,02) denote two 
solutions with the same initial data which exist on some common time interval 
[0,T] and let g = fi — f2, -0 = 0i — '/'2 ■ For brevity we also let X{t) = \\g{t)\\aa + 
\\Dtp{t)\\oo- Using the Vlasov equation we obtain 

Sg-[iS^,)p+ (1 +p^)-'/'-V,^i] ■ Vpg 

= 45^01 + [iS^)p+ (1 V^i/;] • Vp/2 + 4/2 Sip. 

The right hand side of this equation is estimated in modulus by CtM^)- Hence 
integration along characteristics implies 

||5(i)lloo <Ct / X{s)ds. 
Jo 

Taking the difference of the representation formulas for the derivatives of <j)i 
and 02 and estimating as in Lemma [S] we obtain 

mit)\\oo<CT I X{s)ds. 

Jo 

By the previous two estimates A = on [0,T], and uniqueness is established. 
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Assume now that the local solution has been extended to a maximal time 
interval [0,Tmax[, and assume further that A(T„iax) <c». To complete the proof 
of Theorem^we have to show that Tmax = oo. Assume that Tmax < oo. Then due 
to LemmaslHElEl andglthe quantities 7^(0, ll/Wlloo, \\(l>{t)\\oc, ||£'0(t)||oo, and 
||^^</'(i)||oo are uniformly bounded on [0,Tmax[- ForanytoG [OjTmaxl we consider 
the Cauchy problem for the Nordstrom- Vlasov system with f{to), (p{to), dt(j){to) 
prescribed as data at time to] let us ignore for the moment the problem that 
these data do not have the regularity required in our local existence result. By 
the same arguments as in Proposition 21 we conclude that this Cauchy problem 
has a local classical solution in an interval [io,io where S>0 can be chosen 
independently of to- To see the latter observe that all the quantities which 
entered the definition of the local existence time T*, cf. H5.2|l . are now uniformly 
bounded on [0,Tniax[- Hence, if we take to sufficiently close to Tmax, the solution 
can be extended beyond the time Tmax, which is a contradiction to its 
maximality. Thus we have shown that the bound on A implies that the maximal 
solution is global. 

To complete the proof let us deal with the technical catch that 
f {to) , (j){to) , dt(t>{to) do not have the regularity required for our local existence 
result. Instead of using Proposition ^directly we define a new sequence of it- 
erates as follows: For {f^°\4'^°'') we take a global extension of {f,4>)\[o,to] with 
the same regularity as required of the solution and with ||/^°HOIIoo, ||</'^°H*)lloo, 
||£'^'^°-'(t)||oo uniformly bounded in t. Given the (n— 1)'^* iterate we define the 
n**^ iterate exactly as before. In particular, all these iterates coincide with the 
solution on the time interval [0,io], and the term (j>o in H3.1|l . which is the one 
in which the "loss of derivative" problem arises, is the one determined by the 
field data at t = Q. It is now straightforward to repeat all the estimates for these 
new iterates and to show that they converge to a solution on a time interval 
[0,to-|-(5] where S can be bounded away from zero if the quantities entering the 
definition of T* in Proposition 0] are bounded away from infinity, cf. H5.2() . To 
put it shortly: Data obtained by evaluating a classical solution at some time to 
are admissible data for the local existence result. 

6 Proof of Theorem [2] 

The proof of Theorem ^ is not affected by changing the sign of the right hand 
side of H1.5|l and so the Cauchy problem for the "repulsive" Nordstrom- Vlasov 
system has a unique classical solution on a maximal interval of time [0,Tmax[- 
Let us give data which are homogeneous in the ball i3i?(0) and assume that 
Tniax > R- Then the solution remains homogeneous in the cone 

n{R)^{{t,x):\x\<R-t,0<t<R}. 

This property can easily be understood as a consequence of finite propagation 
speed, but for the sake of completeness we give a proof of this claim. We use the 
approximation sequence which we define as in the proof of Theorem^ but with 
the sign of the second term in H5.1|) reversed. The pair 

(/(O)^ 0(0)) is Jiojno 

geneous 
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in Q,{R) by definition. Assume this is true at step n—1. If {t,x) Eil.{R), then 
{s,Xn{s)) wiU be contained in il{R) for all s<t. Hence 

= (s,o)p„, {t,x) G n{R),s < t. 

as 

Integrating this equation we obtain 

P„(s)=pexp((/.("-i)(t,0)-(/.(""i)(s,0)) 

and so P„(O)=pexp(0("-i)(t,O)-(?;.i,"(O)). Hence 

/(")(t,x,p) = /"(O,P„(O))exp(40("-i)(t,O)-40i,"(O)), {t,x)€n{R) 

and so /^"^ is homogeneous on n{R). By the Kirchhoff formula H5.1|l with the 
plus sign between the two terms this is true for t/)'^"^ as well. Since the solution 
(/, (j)) arises as a uniform limit of the iterative sequence on some time interval we 
conclude that it is homogeneous in the truncated cone flriR) — {{t,x) : < i? — 
s,0<s<T} for some 0<T<R. Choose T <R maximal such that the solution 
is homogeneous on flf{R). IfT < R then we could construct the approximation 
sequence with data /(T), (p{T), dt4){T) at time T, which are homogeneous on 
the ball \x\ <R — T, and we would be able to prove homogeneity of the solution 
in the truncated cone ^f_f_s (R) for some d>0 which is a contradiction to the 
maximality of T. 

Hence the solution (/,(/)) is homogeneous on the cone ^r{R) = n{R). If by 
abuse of notation we denote this spatially homogeneous solution on Q{R) by 
f — f{t,p), 4>^4){t) then this pair satisfies the equations 

^=^^, (6.1) 

^,(t)^ I f{t,p)^= (6.2) 
J \/l+p2 

dtf-^p-Vj,f = ^f, (6.3) 

for all [0,i?], where ■ = d/dt. Hence 0>O, and since 0(0) >0, is increasing. 
Moreover 

Integrating this inequality we obtain /i(t) > /i(0)e'^^*-'~'^'-*'-'. Hence by IjG.lfl . 0> 

^(O)e'^-'^(o), and 

dt L 

By our assumptions in Theorem |21 the term in square brackets is non-negative 
at time t = 0. Hence (^(t))^ - 2^(0)e-'^(°)e'^(*) > for all [0,i?]. Setting A = 
2/z(0)e~'^(°\ the latter inequality entails 

_2l(e-*/2)>VA, 
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which implies 

g-0(t)/2<g-0(o)/2„i^^^ te[0,i?]. 

Since by our assumptions in Theorem |21 the right hand side of this inequahty 
is not positive for t = R this is a contradiction. Hence the length of the maxi- 
mal existence interval cannot be larger that R, and the proof of Theorem |21 is 
complete. 

7 The one dimensional case 

In this final section we study the one dimensional system (|1.9|) . p.lO|l . 

in particular, x,pGM. and all integrals are one dimensional. We integrate the 

one dimensional wave equation (|1.9() along its characteristics to obtain 

Mit^x) = ^[{dt(b + d^cj))iO,x + t) + {dt(b-d^(b)iO,x-t)] 

1 /•* 

— - / [ii{s,x + s~t)~ ii{s,x — s + t)]ds, (7.1) 

2 Jo 

5,0(i,x) = ^[{dt(b + d,cl))iO,x + t)-{dt(b-d^(b)iO,x-t)] 



2./0 



[^{s,x + s — t) — ii{s,x — s + t)] ds. (7.2) 



Using these and the resulting formulas for (p itself and its second order derivatives 
instead of the representation formulas in Section^it is straight forward to show 
that initial data of the regularity stated in Theorem|3launch a classical solution 
on some maximal time interval [0,Tmax[, and Tniax = oo if (f) and the momenta 
in the support of / remain bounded, i.e., the analogue of Theorem 1 holds. The 
only point that needs to be checked here is that a bound on 7'(i), ||/(t)||oo, 
||(/)(t)||oo, and ||Z30(i)||oo implies a bound on the second order derivatives of (j), 
i.e., the analogue of Lemma holds. The crucial terms to estimate if one for 
example takes another time derivative in (|7.1II are of the form 

{d.f){s,x + s-t,p)^^ (7.3) 

To get rid of the derivative on / we use the relation 

{d,f){s,x + s-t,p) = j^[dp([{Scf>)p+{l+p^)-'/^d,cl,]f) {s,x + s~t,p) 

+ {Scbf){s,x + s--t,p)--ds[f{s,x + s-t,p)] 



which follows from the Vlasov equation (|1.11|) . Upon substitution into (|7.3|l we 
can integrate by parts with respect to p and s respectively. The resulting terms 
will contain only first order derivatives of (p, no derivatives of /, and kernels 
which are bounded as long as p is bounded away from infinity. 



20 



Assume now that T^ax < oo . Then conservation of energy implies that 4> 
remains bounded on the x-support of /. To see this, note that for i?o>0 
sufHciently large, f{t,x,p) = for |x| > i?o + t , and hence we can change (j) outside 
this forward lightcone without affecting the solution inside. Together with the 
bound on the spatial integral of dxtj) coming from conservation of energy this 
implies a bound on cj) on the a;-support of /. By (|2.7|) this in turn implies a 
bound on / as well; e"^"^/ is constant along characteristics, cf. (|2.6|l . Therefore, 

\\Kt)\U<2\\fit)\u T'^ -^<cHi+rit)) 

with C > independent of t. Combining this with the representation formulas 
(ETTJ and ((7^ implies 

m{t)\\oo <C [I + \n{l+V it))]- 

Using this estimate in (|2.3|l of the characteristic system and integrating we 
obtain the inequality 

V{t)<Po + C f (l + -p(r))ln(l + -p(r))dT 
Jo 

which holds for all t G [0,Tmax[ with some constant C > which docs not depend 
on t. Since by assumption Tmax < oo this implies that V is bounded as well 
which is a contradiction. Hence Tmax^oo. 
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